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Introduction 

The masses of compact objects like white dwarfs, neutron 
stars and black holes are fundamental to astrophysics, but very 
difficult to measure. Close binary systems consisting of a visible 
primary and an invisible compact object are very useful to this 
end as the companion mass can be constrained from the radial 
velocity and the light curve of the primary if the primary mass 
and the orbital inclination are known. The latter can be measured 
in systems that are eclipsing or show tidally locked rotation. If 
the mass of the primary can be estimated, the companion mass 
can then be derived. The assumption of tidally locked rotation 
has often been used to determine the masses of neutron stars and 
black holes in known X-ray binaries (see Charles & Coe 120031 
for a review). The same technique has recently been applied to 
KPD 1930+2752, a short period binary consistingof a sublumi- 
nous B stars and a white dwarf (Geier et al. l2007l H 

sdB stars are core helium burning stars with very thin hydro- 
gen envelopes and masses around O.5M (Heber et al. 1 1 9861 ). A 
large fraction of the sdB stars are members of short period bi- 
naries (Maxted et. al 1200 It Napiwotzki et al. 120041 ). For these 
systems common envelope ejection is the most probable forma- 
tion channel (Han et al. 2003). The companions of sdBs in these 
systems are predominantly white dwarfs implying that the sys- 
tem has undergone two phases of common envelope ejection. 

In order to establish bound rotation for an sdB star the syn- 
chronization time scale has to be smaller than its evolutionary 
life time (?ehb ~ 10 8 yrs). Two theoretical concepts to com- 
pute synchronization times have been developed by Zahn ( 1 19771 ) 
and Tassoul & Tassoul (1992), respectively. While the mecha- 
nism proposed by Zahn is not efficient enough to fully fit the 
observed levels of synchronization, the more efficient one of 
Tassoul & Tassoul is matter of a controversy, given its free pa- 
rameter dependence (see Claret et al. 1995, 1997 and references 
therein). Especially in the case of hot stars with radiative en- 
velopes, where tidal forces are less effective in synchronizing the 
stars, the predictions of the two theoretical models at hand can 
differ by orders of magnitude. All prior studies were undertaken 
to match observations of hot main sequence stars. Hot subdwarf 
stars have similar temperatures as B-type main sequence stars, 
but are much smaller and the internal structure of these helium 
core burning objects is different. In addition, the fraction of sdBs 
residing in close binary systems is among the highest known of 
all types of stars. Observations of hot subdwarfs could provide 
a new benchmark to study the yet unresolved problem of tidal 
dissipation in radiative stellar envelopes. 

Independent observational constraints are needed to prove or 
disprove synchronized rotation in hot subdwarf stars. Ellipsoidal 
variations can be used to verify synchronization of the stellar 
surface because the light variations would then have to occur at 
exactly half the orbital period. Two sdB + white dwarf binaries 
are known to show ellipsoidal variations at half of the orbital pe- 
riod (KPD 0422+ 5421 Orosz & Wade |T999l KPD 1930+2752, 
Geier et al. l2007b . However, both systems have short orbital pe- 
riods of about 0. 1 d and theory predicts synchronization times 
much smaller the evolutionary time scale. As the synchroniza- 
tion time strongly increases with increasing period, we expect 
an upper limit to the period to exist at which the assumption of 
tidally locked rotation breaks down. To this end it would be of 
utmost importance to find ellipsoidal variations in an sdB binary 



of longer period and to provide a stringent test for the theory of 
synchronization. 

Recently a suitable object has been found. PG 0101+039, 
an sdB+WD binary (P=0.567 d, Maxted et al. I200TT ) was dis- 
covered to show very weak luminosity variation at half the or- 
bital period in a light curve in a 16.9 day long, almost uninter- 
rupted light curve obtained with the MOST satellite (Randall et 
al. l2005l 

In order to verify that we indeed see ellipsoidal variations, 
we have to show that the observed light curve can be consis- 
tently modelled. Beforehand, we have to derive the complete set 
of system parameters. As the spectrum is single lined, the anal- 
ysis of the radial velocity curve yields the mass function only. 
Complementing it with an estimate of the sdB mass and with 
measurements of the sdB's projected rotational velocity as well 
as its gravity allows to solve for all binary parameters and com- 
pute the light curve. 

1. Binary parameters 

1. 1. Radial velocity curve 



Based on spectra obtained in 1998, Moran et al. ([1999) deter- 
mined the period P = 0.569908 ± 0.000007 d. However, these 
ephemerides are not accurate enough to phase the MOST pho- 
tometry because the time span of six years between spectro- 
scopic and photometric observations is too long. Therefore we 
combined the velocities of Moran (1999 ) with those from eight 
MMT -spectra taken in 1996, 1997 and 2002 (Randall et al. 
2005 ) and five spectra obtained with the Steward 2.3m Bok tele- 
scope from 2000. The latter were determined using the double- 
precision version of the IRAF fxcor package, against the com- 
bined template for the star. In addition we obtained three high 
resolution spectra (R = 30000) in 2000 at the 2.2 m telescope 
at the German-Spanish Astronomical Center (DSAZ) equipped 
with the FOCES-spectrograph. RVs were determined by cross 
correlation with a model spectrum at rest wavelength. 

This provided us with 57 velocities distributed between 1996 
and 2000. A sine curve was fitted to the observed velocities us- 
ing an x 1 minimising method and the power spectrum was gen- 
erated (Napiwotzki et al. 1200 IK The orbital parameters were 
measured: y = 7.3 + 0.2 km s" 1 , K = 104.5 ± 0.3 km s" 1 , 
P = 0.569899 + 0.000001 d and HJD(0) = 2452545.064275. 
The period error was derived with a bootstrapping algorithm. 
The period is sufficiently accurate to allow the MOST photome- 
try to be phased properly. However, it should be kept in mind that 
the error quoted is a purely statistical one. To investigate pos- 
sible systematic errors, which can occur when radial velocities 
obtained with different instruments are combined, we also used 
another method to verify our results. In this case the fitting was 
performed with the system velocity y as additional free parame- 
ter for every individual dataset. The results using both methods 
were perfectly consistent within the quoted errors. However, ad- 
ditional systematic effects may still be present. 

1.2. Gravity and projected rotational velocity 



1 The companion is so massive that the system mass may exceed the 
Chandrasekhar mass, making the system a viable Supernova la progen- 
itor candidate in the double degenerate scenario. 



Low resolution spectra obtained by Randall et al. (2005J and 
Maxted et al. (2001) were used to derived the atmospheric pa- 
rameters. Particular attention should be paid to the gravity deter- 
mination as it provides a mass-radius-relation and its error prop- 
agates into the mass determination (see. Sect. 4). 

Synthetic line profiles calculated from metal line-blanketed 
LTE model atmospheres with solar metal content (Heber et al. 
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Fig. 1. Companion mass as a function of primary (sdB) mass. 
The dotted vertical lines mark the theoretical sdB mass range for 
the common envelope ejection channel (Han et al. 2003 ). The 
dashed lines mark the error limits of the companion mass. 



2000) were matched to the observed Balmer and helium line 
profiles using a x 1 fit procedure described by Napiwotzki et al. 
( I1999I ). The resulting parameters are r er r = 27 700 K, log g = 
5.55, logN(Ue)/N(H) = -2.62 from the Randall et al. spectra 
and r eff = 27 300 K, log g = 5.50, log N(Ue)/N(U) = -2.71 
from the Maxted et al. spectra with formal statistical fitting er- 
rors of less than 100 K, 0.02 dex and 0.02 dex, respectively, 
which are unrealistically low. The true uncertainties are domi- 
nated by systematic inaccuracies in both the observations and 
model atmospheres and can be estimated from repeated obser- 
vations and the use of different model grids. Taking into ac- 
count the discussion of typical systematic errors applying this 
method in Geier et al. (2007 ) we adopt r eff = 27 500 ± 500 K, 
logg = 5.53 + 0.07, logN(He)/N(U) = -2.66 + 0.1. 

In order to derive v rot sin i and the elemental abundances, we 
compared the observed high resolution spectra with rotationally 
broadened, synthetic line profiles. The projected rotational ve- 
locity was measured simultaneously with the elemental abun- 
dances to v rot sin « = 10.9 + 1.1 kms 1 using 17 suitable metal 
lines. 



1.3. Analysis 

The analysis strategy is the same as for KPD 1930+2752 and 
therefore is described here only briefly. For details we refer the 
reader to Geier et al. (2007*1). 

Since the spectrum of PG 0101+039 is single-lined, it con- 
tains no information about the orbital motion of the companion, 

and thus only the mass function / m = ,.^" m '' s ' n '.. 



= can be 



calculated. Although the RV semi-amplitude K and the period P 
are determined by the RV curve, M s( jb, M comp and sin 3 i remain 
free parameters. 

Nevertheless, the masses can be constrained by assuming 
tidal synchronization. Combining the orbital parameters with 
an estimate of the sdB mass and with the measurements of its 
v rot sin i and gravity allows the mass of the invisible companion 
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Fig. 2. Lightcurve data with superimposed models. The two solid 
curves confine the best fit models with parameters and associ- 
ated uncertainties derived from spectroscopy under the assump- 
tion of orbital synchronization and M s[ jb = 0.3 M . The two dot- 
ted curves confine models for M sr jB = 0.7 M . Phasing the data 
to the orbital solution results in a significant phase shift (upper 
panel). In the lower panel the model lightcurve is shifted by -0. 1 
in phase. 



to be constrained tightly. The mass of the sdB primary is con- 
strained from the population synthesis models (Han et al. 2003 ) 
which predict a mass range of M sc ib = 0.30-0.48 M for sdBs 
in binaries, that experienced a common envelope ejection. The 
mass distribution shows a sharp peak at a mass of about 0.47 M . 
If the companion is synchronized the rotational velocity 
can be calculated. The radius of the primary is given 



by the mass radius relation R = yj M ^ G ■ The measurement of 

the projected rotational velocity v 10t sin i therefore allows us to 
constrain the inclination angle i. With M s( jb as free parameter the 
mass function can be solved and the inclination angle as well as 
the companion mass can be derived. Because of sin i < 1 a lower 



limit for the sdB mass is given by M s( jb > 



4tt 2 G 



There are no 



spectral signatures of the companion visible. A main sequence 
companion with a mass higher than 0.45 M can therefore be ex- 
cluded because its luminosity would be sufficiently high to be 
detectable in the spectra. 

The relation between the primary and the secondary mass 
is shown in Fig. Q] The allowed mass range for the companion 
is consistent with that of a white dwarf and therefore consistent 
with the common envelope ejection scenario. For the most likely 
sdB mass of 0.47 M the binary parameters are: R s ^b - 0.19 + 
0.02 fl , M WD = 0.72 + 0.10 M , inclination i = 40 + 6° and 
separation a = 3.1 ± 0.4 R Q . 
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2. Light curve and elipsoidal variations 

The MOST photometric data were folded on the orbital period 
and re-binned to get a better signal to noise ratio (see Fig. [2] As 
PG 1 1 +039 is a pulsating sdB star of V 1 093 Her type the light 
curve had to be pre-whitened for the pulsational frequencies be- 
forehand. 

Each bin contains more than 400 original measurements. 
Light curve variations at half the orbital period with a semi- 
amplitude of 0.025% were detected. Therefore PG 0101+039 
shows the smallest ellipsoidal variation ever measured. In or- 
der to compare with synthetic light curves we had to properly 
phase the photometry taking the zero point from the orbital solu- 
tion. Because the amplitude of the variation is very low, Doppler 
boosting arising from the orbital motion affects the lightcurve 
significantly. A resulting factor of (1 - v(t)/c) was applied to the 
total flux to correct for this effect. 

The light curve was modelled with the light curve synthe- 
sis and solution code MORO based on the model by Wilson 
& Devinney (|1971| l. The details of the Bamberg implementa- 
tion are given by Drechsel et al. ( 119951 1. The software uses a 
modified Roche model for light curve synthesis. It is capable of 
simulating the distortions of the stars induced by a companion. 
Light curves for different component masses and orbital inclina- 
tions were synthesised. For typical values of M s dB appropriate 
ranges of sin i and values were computed as described in 
Sect. [TJcovering the full parameter space (including error limits) 
allowed by the spectroscopic analysis. 

We used sdB masses from 0.3 - 0.7 M Q for the model light 
curves and compared it to the observations in Fig. [2] The syn- 
thetic light curve matches the semi-amplitude of the observed 
one quite well. Taking into account the extremely low ampli- 
tude of the variation, the consistency with the model is remark- 
able. However, there is a significant phase shift between the ob- 
served and the predicted light curve. The best fit to the data is 
phase shifted by -0.1 with respect to the model calculated with 
proper orbital phase (see Fig. |2). Despite our efforts to derive 
a high precision orbital period, we can not rule out completely 
that this shift is caused by systematic effects when RV measure- 
ments from different instruments are combined. To match the ob- 
served light curve would require the period to differ by 0.00008 d 
from our results. Given the overall consistency of our orbital pa- 
rameter determination, such a large deviation (80x period error) 
seems to be rather unlikely. 

3. Discussion 

Tidally locked rotation in close binary systems has been assumed 
to measure masses of invisible compact companions, in particu- 
lar in X-ray binaries. The synchronization time is very difficult 
to calculate for stars with radiative envelopes and plagued with 
large uncertainties. Therefore, observational constraints are of 
utmost importance. Ellipsoidal variations can be used to verify 
the assumption at least for the surface layers. We applied this 
technique to the sdB/WD binary PG 0101+039, for which a very 
weak luminosity variation at half the orbital period has been dis- 
covered in a light curve in a 16.9 day long, almost uninterrupted 
light curve obtained with the MOST satellite. 

From spectroscopy we measured the mass function, appar- 
ent rotation and surface gravity of PG 0101+039. Stellar evo- 
lution models suggest that the sdB mass is close to 0.47 M . 
Assuming tidally locked rotation, this information is sufficient 
to solve for all parameters of the binary system. The companion 
mass is found to be Mwd = 0.72 ± 0.10 M , typical for a white 



dwarf. The light curve was then synthesised and was found to 
match the observed amplitude well. However, a problem with 
the phasing of the light curve to the radial velocity curve became 
apparent. Due to a six year difference between the MOST pho- 
tometry and published radial velocities, the phase errors were far 
too large for any conclusion to be drawn. Therefore we added 
16 radial velocities from three observatories. The statistical er- 
ror of the period decreased sufficiently to enable proper phasing 
of the photometry. The synthesised light curve was found to be 
offset by 0.1 cycles from the observed one indicating that our 
systematic error estimate may be overly optimistic. Alternative 
explanations like supersynchronous rotation of the sdB that may 
cause the observed phase shift seem to be unlikely because a de- 
viation of 10% from equilibrium would require fast rotation of 
the sdB. In this case the inclination would be very low and the 
companion mass would rise dramatically. 

A simultaneous measurement of the radial velocity curve and 
a high precision light curve would be necessary to solve this 
problem since the theoretical understanding of angular momen- 
tum transfer in hot stars with radiative envelopes is still very 
limited. In conclusion, we found strong indication that the sur- 
face rotation of the sdB star PG 0101+039 is tidally locked to its 
orbit. 

The synchronization times for any given type of primary 
depend strongly on the orbital period (Zahn 119771 Tassoul & 
Tassoul 1992). Hence, other sdB stars in close binaries should 
also be synchronized if their orbital period is less than that 
of PG 0101+039 (P = 0.567 d). Hence we conclude that tidally 
locked surface rotation is established in sdB binaries with orbital 
periods of less than half a day. Hence the assumption of tidally 
locked rotation can be safely applied to such systems, even if 
they are too faint to measure such extremely small light varia- 
tions as observed here. 
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